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Abstract 

We present a field theory of Jain's composite fermion model ffl, as generalised 
to the bilayer quantum Hall systems. We define operators which create com- 
posite fermions and write the Hamiltonian exactly in terms of these operators. 
This is seen to be a complexified version of the familiar Chern Simons theory. 
In the mean-field approximation, the composite fermions feel a modified ef- 
fective magnetic field exactly as happens in usual Chern Simons theories, and 
plateaus are predicted at the same values of filling factors as Lopez and Frad- 
kin [|[and Halperin |J. But unlike normal Chern Simons theories, we obtain 
all features of the first-quantised wavefunctions including its phase, modulus 
and correct gaussian factors at the mean field level. The familiar Jain rela- 
tions for monolayers and the Halperin wavefunction for bilayers come out as 
special cases. 
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I. INTRODUCTION 

With the development of techniques for growing GaAs heterostructures containing two 
separated layers of two-dimensional electron gas, experimental work on the quantum Hall 
effect has been extended to such bilayer systems as well (see for example references |l| . 
01) @) 01 )• New plateaus in Hall conductivity have been observed at filling fractions not 
seen in single layers such as at v — |. On the theoretical front, a large body of work has 
already been done on bilayer systems. An extensive list of references to this literature has 
been given in the lucid review of this subject by Girvin and MacDonald |J and in the paper 
by Moon et al ||. In particular, as background and motivation for our present work let us 
recapitulate the following theoretical developments : 

(i) A major step which helped in the study of quantum Hall effect was the proposal 
of elegant and yet very accurate first-quantised N-particle wavefunctions for the ground 
state and quasiparticle excitations at the Hall conductivtity plateaus. This was pioneered 
by Laughlin for the mono layer case with his famous wavefunction for the filling fractions 
°f v = 2m+i E3' F° r double layer systems (viewed as a two component system carrying 
a pseudo-spin layer index ) a generalisation of the Laughlin wavefunction was proposed 
long ago by Halperin || . The Halperin wavefunction ip mitm2>n is labelled by three integers 
mi > m 2 an d n — of which m\ and m 2 must be odd - which determine the filling fraction v. 



For example, it was proposed by Yashioka, McDonald and Girvin [jTTJ that the plateau at 
v =\ seen in the bilayer system corresponded to the wavefunction ^3,3,1- 

(ii) For mono-layers , the fractional quantum Hall effect plateaus and their phenomeno- 
logically very successful wavefunctions were derived or justified from underlying composite 
particle formation postulates. Jain |l| presented a theory of flux-electron composite fermion 
formation. Jain's theory related Hall plateaus and their electron wavefunctions at fractional 
fillings to corrresponding plateaus and wavefunctions of the composite fermions at integral 
fillings. 

(m) On a different front, the fractional effect was studied in field theoretic formulations. 



Based on the observations of Girvin and MacDonald [[L2| that these systems seem to exhibit 
off-diagonal long range order, Zhang, Hansson and Kivelson [|H| constructed a Chern Simons 
field theory of a Landau Ginsberg order parameter field for quantum Hall effect at fractional 
fillings v = 2m+i ■ Their bosonic order parameter field corresponds to composites of electrons 
with an odd number of fluxons. A similar Chern Simons field theory, but with even -integer 



coupling, was studied by Lopez and Fradkin [|14| . This corresponds to having composites 



of even number of fluxons with electrons and gives a field theoretic formulation of Jain's 
theory. Subsequently, Lopez and Fradkin |2| also extended their fermionic Chern Simons 
field theories to the bilayer case and predicted possible Hall plateaus for a large family of 
filling fractions in each layer. For related work on partially polarised electrons see Mandal 



and Ravishankar ||15|. Ezawa and Iwasaki [16| also studied bilayer systems through a bosonic 



Chern Simons theory. They solved it using self-duality equations which hold in the absence 
of Coulomb interactions, which they treat perturbatively in the short distance limit. 

(iii) All these papers in references [[13[] , |IB| , |1"4|| and |2| extracted first quantised N- 
particle wavefunctions from their field theoretic ground states. But in the mean field ap- 
proximation, only some aspects of the Laughlin wavefunction emerged. In the Zhang et al 
theory [DJ, the Landau Ginsberg field incorporates only the correct phases of the electronic 



correlations and not their modulus. The all important zeroes in the Jastrow correlation 
factors as well as the gaussian factors in the Laughlin wave functions emerge only upon 



including fluctations about the mean field theory. In ref. [jIJJ the modulus of the Laughlin 
wavefunction is derived by a very different and ingenious method starting from the low-q 2 
limit of correlation functions. But such a result holds only in the long distance limit. The 
same is also true of the Lopez -Fradkin bilayer work where again they have obtained 
the Halperin-Jain wavefunction for some cases, but only in the long distance limit. The 
wavefunctions obtained in the Ezawa-Iwasaki work hold only at short distances and are in 
the absence of Coulomb interactions. 

In this paper we present a modified Chern Simons field theory for bilayer systems, which 
yields naturally at the mean field level Jain's model of composite fermion formation as 



generalised to bilayer systems as well as the corresponding wavef unctions. (Jain's well known 
monolayer results also come out as a special case.) The present work is an adaptation 
to bilayer fermionic composite operators of our earlier work with Sondhi [[Tj]] where we 
had presented an exact field theory of the Read operator |TSfl . There we had employed a 
complexified version of the Chern Simons Landau Ginsberg field theory which enabled us to 
reproduce all the features of the Laughlin wavefunction already at the mean field level. We 
use an appropriate generalisation of the same method here. We explicitly construct operators 
that create bilayer composite fermions, using a non-unitary transformation acting on the 
parent electron field operator . Exact anti-commutation rules and an exact Hamiltonian 
are written in terms of these composite fermion operators. In the mean field approximation 
this Hamiltonian relates electrons at fractional fillings to composite fermions at integer 
fillings. These fractions and integers are seen to be related precisely by the formulae given by 
Lopez and Fradkin. Equations are also obtained akin to Jain's , but generalised to bilayers, 
which relate electron wavefunctions to composite fermion wavef unctions. The Halperin 
wavefunctions are obtained as specific examples. 

In our theory, all aspects of these wavefunctions - the phases and moduli of their Jastrow 
correlations, and appropriate gaussian factors - emerge in the lowest order of the mean-field 
approximation. No short distance or long distance approximation is used nor any lowest 
Landau level restriction put in by hand. That our method of reference |l7l can be used to 



generate composite fermions was also pointed out by Wu and Yu [19] for the monolayer case. 



II. COMPOSITE FERMION OPERATORS 

Consider a double layer of two-dimensional electrons of mass /x and charge e, placed in a 
uniform perpendicular magnetic field of strength B which corresponds to a vector potential 
in the symmetric gauge of A(x) = \B k x x, where k is the unit vector perpendicular to the 
plane. We take the electrons to be fully spin polarised along the B field and hence suppress 
spin for simplicity. Suppose their interaction potential is V(x — x') and the scalar potential 



Aq represents any uniform and/or random impurity electric fields in the problem. Let ^ a (x) 
denote the bilayer electron quantum field with a = 1,2 standing for the layer index. It 
obeys the equal-time anticommutation relations, 

{*«(*), Ma?)} = {*£(*), ^W )} = ° 

{*«(£), ttj(f)} = 5 a>p 5^(x-x') . (2.1) 

Clearly ty a is the full electron field and not just its lowest Landau level projected part. Note 
that unless explicitly specified, it is understood throughout this article that repeated indices 
a, (3 etc. are not to be summed over. The second quantized Hamiltonian that describes our 
system is, 

H = Jd 2 xJ2 ^i(x)(^—D 2 + eAo)^ a (x) + ]T - J J d 2 xd 2 x 5p a (x) V af3 (x - f ) 5 Pf3 (x 



a,(3 

(2.2) 



Here, D = V — i-£-A is the covariant derivative inclusive of the vector potential of the 
uniform magnetic field B . p a {x) = ^' a (x)^ a (x) is the electron density operator in the 
layer a whose deviation from its mean value p a is Sp a (x). 

Next we define our operators for the composite fermion field Xa{x) and its canonical 
conjugate field U a (x), by 

Xa(x) = e~ M£) ^ a (x) 

U a (x) = *i(5)e J «W, (2 .3) 



where 
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J a (x) = Y^R a p J d 2 x' [ Pl3 (x') \og(z - z')] - rj-^-, (2.4) 

In the above equation r a and the 2x2 matrix R a p consist of some numbers which will be 
fixed shortly, z = x\ + 1x2 is the complex coordinate on the plane and I = J ^ is the 
magnetic length. Notice that Xa and n a are not hermitian conjugates of each other since 



J a has both hermitian and anti-hermitian pieces. Nevertheless, as we now show, they form 
a pair of canonically conjugate Fermi fields, one in each layer. 

Note that the only operator appearing in J a (x) is the electron density pp(x) = 



typ(x)typ(x), which obeys the commutation relation, 



\p a (x),V(i(#)\ = -8 a p V a (g) 5 2 (x-x') (2.5) 

Therefore the following identities follow : 

e --M*0 m p (x') = (z- z') Ra0 *p{&) e~ Jo(a) 
^l(x') e~ J ^ ] = (z- z ') Ra0 e- J <*W ^l(x') . (2.6) 

Using these identities one can verify that 

Xa (x)xp(x') = (-l) R ^ +1 (z - z') Ra0 ~ R " a Xp(x')Xa(x') (2.7) 



We can see that the explicit z-z' dependence in ( |2.7|) drops out if R a p is chosen to be a 
symmetric matrix. Further, if its diagonal elements are even intergers (Rn = 2s\ and 
i?22 = 2S2) then the field x i n each layer anticommutes with itself , as desired of composite 
fermions. The off diagonal element Ru = R21 can be taken to be an integer n. Depending 
on whether n is odd (even), the fields at two different layers will commute (anticommute). 
In short the requirement that the composite fields defined in eq ( |2.3|) be fermi fields restricts 
the matrix R a p to have the form 

/2si n \ 
R a p = (2.8) 

V n 2s 2 ) 

exactly in accordance with ref @. 

The same choice of R a p also yields the canonical anti-commutator between Xa and IL3. 

We have, upon using the identities (|2.6| ), 



X a {x)Tlp{x') = e-^^f^fV^ 

= (_i)«^+i Up(x') Xa (S) + 8 a p5 2 (x-x') . (2.9) 



Thus, since R\\ and R22 have been even integers, the fields \a and Ii a form a pair of mutually 
conjugate Fermi fields for each layer a. This is despite the presence of non-unitary factors 
in their definition in Eq. ( |2.3| ). However, in contrast to standard fermi field theories, here 
Ha is not equal to xL- Instead they obey the more complicated relation 

n Q (x) = xUx)e MS)+jta{ * ] . (2.10) 

This has to be borne in mind in doing manipulations with these composite fermion fields. 
We will define the composite fermion density p a by 

p a (x) = U a (x)x a (x) . (2.11) 

The corresponding number operator N a = J d 2 x p a satisfies 

[N a ,U /3 (x)} = S af3 U a (x), (2.12) 

i.e. the operator H a (x) creats one extra composite fermion in the a th layer. Notice that we 
have used the same symbol p a (x) for the this composite fermion density as we did for the 
original electron density since the definition given in ( |2.11| ) also satisfies 



p a {x) = ¥ a {x) * a (x) (2.13) 

It should be emphasized that our composite fermion operators Tl a (x) and Xaip) are 
defined over the same space-time domain as the original electron field, as is natural in any 
field operator transformation in a field theory. This means that the area of the Hall sample 
is the same, whether we consider electrons or composite fermions. Hence, since the densities 
of both types of fermions has been shown to be the same, the total number of composite 
fermions is also the same as the number of the original fermions. Note that in obtaining the 
identities ( |2.6| ) we have used the expression ( |2.13| ) for the density in terms of the electron 
field and the associated commutator. If we were instead to use the expression ( |2.11| ) in 
terms the composite fermion operators and the associated commutators, we can analogously 
obtain the identities 



e- J "W xp tf) 



[z-z') Ra& Xp{&) e ~ M£) 



U p {x') e~ MS) = {z- z ') R ^e- MS) U (x' ) . 



(2.14) 



Having denned the composite fermion fields and obtained their commutation relations, 
let us next rewrite the Hamiltonian fl2.2| ) in terms of them. First consider the covariant 
derivative on the electron field. We have, 



DV a (x) 



D(e J ^ Xa (x)) 

(V - i^-/{x)) (e Jaig) Xa (x)) 



?Ja{x) 



(V - i—A(x) + VJ a (f)) Xa(x) 



he 



Ja{x) 



IC 



(D ~ T-V<*{x)) Xa(x) 



(2.15) 



where, 



v a (x) = i — VJ a (x) . 



(2.16) 



Hence, 



D^ a = e J « (D - 1—VaYXa 

he 



(2-17) 



Inserting this into the starting Hamiltonian ( |2.2j ), and using Eqs 

-h 2 



and ( |2.11|) we get 



H 



d x 



UJx) 



le 



■>i ^"¥c ( A + v ^ 2 + eA ° ) Xa(f) 



1 
+ 2 



d 2 x d 2 x' Sp a (x) V a/ 3(x — x') Spp^x') 



(2.18) 



This Hamiltonian in terms of the composite fermion fields defined in ( |2.3|) is exactly equal 
to that of our original electron problem . No approximations have been made so far. Clearly 
this is very similar to a Chern Simons theory but it is more than just a direct generalisation 
to bilayers. As in normal Chern Simons theories, the vector field v a appearing in ( |2.18|) 
above is also constrained in terms of the density by Eq. ( |2.16| ), where J a (x) is defined in 
Q2.4p. But since this J a (x) involves more than just the phase of [z — z'), the field v is not 
the bilayer statistical Chern-Simon gauge field used, for instance, in @. Because J a (x) has 
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real parts, v a is a complex vector field. However, v a will turn out to be simply related to 
Chern Simons fields. 

Let us define a Chern-Simons field for each layer index a by 

— Tic -* f 
S a (x) = V x / <Px? p a {tf)lmlog(z - z') , (2.19) 

or equivalently 

b a (x) = V x a a {x) = - (pop a {x) (2.20) 



where 0o = — is the flux quantum. Following reference JL7), use the Cauchy-Riemann 
condition 

V(Relog2) = V(Imlogz) x k (2.21) 

where k is a unit vector perpendicular to the plane. Using this we get, 

v a {x) = — vJ a {x) 
e 

iTic -* ( t I z I \ 

= — V * { Y, R <*P J d2x ' [pp(x')(Re\og(z - z') + i\m\og(z - z'))\ - ^^f} 

= Y, R ap\ap(x) + i kx ap(x)J -i^ a —— . (2.22) 



III. MEAN FIELD APPROXIMATION 

Thus far everything is exact. Let us now introduce the mean field (MF) approximation, 
by replacing the actual space dependent density operator pp{x) in (|2.22|) by its average value 
pp = ^r— where vp is the filling factor in the layer (3. Under this approximation eq.( [2.20| ) 
becomes 



'aJMF — V X {Cl a ) MF 
= - (poPa 

= -u a B (3.1) 

9 



or equivalently, 



\fla)MF — V a A 



(3.2) 



where A is the applied external vector potential. Then, (|2.22|) reduces to 



- A ^RapVp + i (k X A)( -J2 R »f3Vf3 + r a ) 



(3 * a 

Here we have used the fact that in our symmetric gauge 



(3.3) 



x 



Ant 1 



k x A 



(3.4) 



Then the covariant derivative in the composite fermion Hamiltonian fl2.18|) becomes, upon 



using the MF approximation ( [3. 3D , 



D r 






te 



= V- — All - £iV/>) + i(kxA)l-J2R a0 u + r a ) (3.5) 

Recall that the numbers r a were introduced in our definition of the composite field operators 
r3D and (|2.4|). Thus far we had left them unspecified, but let us now choose them to satisfy 



r a = ^R a pi>p 





(3.6) 



Then the covariant derivative simplifies to 



D n 



ze 

V- — A[ 1 

nc 







(3.7) 



Therefore in the MF approximation, our composite fermion field experiences a Hamiltonian 

-h 2 



H 



d x 



n a (f 

i 



+- 



d 2 X d 2 x' 5 p a (x) V a /3 (x — x') 5pp (x') 



(3.8) 



where 



A 



1 - X! ^^ A 



(3.9) 
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B* a = curlA a = 1 - £i? Q/3 ^ B (3.10) 

V P J 

is the effective magnetic field felt by the composite fermions of the ath layer. Since the 
filling factor is inversely proportional to the magnetic field, this reduced effective magnetic 
field amounts to a correpondingly enhanced filling factor for the composite fermions given 
by 



-1 ~ ' Up RafiVp, 



When these equations are inverted we get 



AV pi 

where 

1 \ / 1 



1 ( 1 

AV p 2 

v 2 = — (n 2si) (3.12) 



A = n 2 - ( — + 2si)( — + 2s 2 ) (3.13) 

\Pl / Vp 2 / 

These are precisely the filling factors obtained for the bilayer system by Lopez and Fradkin 



IV. FIRST-QUANTISED WAVEFUNCTIONS 

The whole purpose of our defining composite fermion operators is to give a field theoretic 
version of Jain's theory of fractional quantum Hall effect as generalised to bilayer systems. 
Following Jain's philosophy, if the effective filling factors p a are integers then one may intiu- 
tively expect the composite fermions to form incompressible quantum Hall ground states. 
In terms of electron coordinates the same quantum Hall state would appear at fractional 



filling factor u a related to the p a by equations (|3.12|) . 

Further, since the composite fermion field Xa can be expressed in terms of the original 
electron field \1/ Q through equation (|2.3|) , if one knows the first quantised wavefunction of 
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composite fermions in a given state, the corresponding wavefunction of electrons in the same 
state is directly obtainable. 

Let (f> Pl , P2 (zi, z 2 , , Zn x ] wi,W2, , wn 2 ) represent the composite fermion wavefunc- 
tion at composite fermion filling fractions p a = ^r, where B* is the effective magnetic 
field felt by the composite fermions in layer a as given in in eq ( |3.10| ). The Z{ are complex 
coordinates of the composite fermions on the first layer and the Wi the complex coordinates 
in the second layer. Note that iViand N 2 , which stand for the number of composite particles 
in each of the two layers respectively, are also the number of electrons in each of the two 
layers, since the density and Number operators of the electrons and those of the composite 
fermions are equal to one another (see eq |2.13| ). We recall that both the electron and the 



composite fermion operators are defined in the same two-dimensional space with the same 
area (sample size). This first quantised wavefunction can be written in terms of the field 
theoretic states and field operators Xa as follows : 

0pi,p 2 (^i)^2, ,z Nl ; w 1 ,w 2 , ,w N2 ) 

= (0|xi(«i)Xi(^)-Xi(«Wi) X2(u)i)X2(w 2 ) X2(w N2 )\MF) (4.1) 

where \MF) stands for the mean-field ground state of the composite fermion system at filling 
factor p a and (0\ is the vacuum state. Meanwhile corresponding to the same state \MF) 
the first quantised wavefunction of the electrons (whose filling fractions in the two layers are 
respectively v^Vz) is given by 

ipv uV3 (zx, Z2, ,z Nl ; w 1 ,w 2 , ,w N2 ) 

= {0\* 1 (z 1 )Vi(z 2 )...*i(z Nl ) ^ 2 ( Wl )^ 2 (w 2 ) V 2 (w N2 )\MF) (4.2) 

To relate the two wavefunctions, one only needs to write the operator if> a in terms of Xa 
using ( P75| ) . We get 

ip„ ljV2 (zi, z 2 , , z Nl ; Wi,w 2 , ,w N2 ) 

= {0\e J ^xM)e J ^Xi{z2)...e J ^ z ^ ) Xi{z Nl ) 

e J ^ X2 (w 1 )e J ^ X 2(w 2 ) e h{wN ^x 2 {w N2 )\MF) (4.3) 
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Next bring all the e J factors in the above expression to the left by commuting them 
across the operators x using the commutators ( [2.14| ). We get 



ipvi,v 2 (zi, z 2 , ,z Nl ; Wl,W 2 , ,w N2 ) 

JVi N 2 N-l N 2 

= R(Zi - Zj) Rl1 I1K - Ml)* 22 UUiZr ~ W S 



,-Rl2 



i>j k>l r s 

Ni N 2 

( 0\ex P {Y.Uzi) + E^K)) 
1 1 

X l(*l)Xl(*2)-..Xl(ZiVi) X2(^l)X2(^2) X2(w N2 )\MF) (4.4) 

Next apply the operators e J to the left on the vacuum state. Notice from its definition in 
eq. (|2.4|) that the only nontrivial operator contained in J a is the density, in the first term. 
The second term in J is just a c- number. In an interacting field theory, the vacuum is not 
generally an eigenstate of the density operator. However in the spirit of the mean field 
approximation being employed in this section, one can replace the density operator by its 
mean value. The mean density of the vacuum is zero. Thus when the operators e Ja in the 
above equation act on the left on the vacuum state the density dependent first term of J a 
can be taken as zero and only the second ( c-number ) survives , giving gaussian factors. 
Therefore 



Ipvuviizii z 2, ,z Nl ; w 1 ,w 2 , ,w N2 ) 

Ni N 2 N x N 2 

= HiZi - Z,) 2S1 l[(w k - W^HHiZr ~ W S 

i>j k>l r s 

_l JVi N 2 



jnin EN 2 + ^ EN 2 ) 

m 1 l 



exp 
( 0|Xi(3i)Xi(3!)-Xi(3Ni) X2(w 1 )x2(w 2 ) X2(w N . 2 )\MF) (4.5) 



In the above equation we have also inserted the matrix elements of R a p from eq (|2.8|) . In 
terms of the composite fermion wavefunction defined in eq. (|4.1| ) we thus get 



VW 2 0*1,^2, ,ZNi ; Wi,W2, ,W N . 2 ) 

Ni N 2 JVi N 2 

= n( Zi - Zj ) 2s > n(w k - Wl y s * nn(^ - «o b 

i>j k>l r s 
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exp 



JVi 



N 2 



jn(nJ2\ z i\ 2 + ^E 

<±l 1 1 



»'; 



Pl,P2\Z~l-l ^2) 



2ATi! WuW 2 , 



w N2 



(4.6) 



This is just the generalisation to double layers of Jain's formula relating wavefunctions of 
electrons at certain fractional fillings to corresponding wavefunctions of composite fermions 
at other related fillings. 

V. DISCUSSION 



Several features of the result ( |4.6| ) are worth pointing out. 

(a) Although this relation holds for any electronic filling u a and the corresponding com- 
posite fermion filling p a , Jain's theory pertains to cases where the p a are integers. Then, 
in the non- interacting limit, the composite fermions will completely fill an integer number 
of Landau levels, giving rise to an energy gap. Therefore the usual arguments in the quan- 
tum Hall literature can be invoked to expect that even in the presence of e-e interactions, 
impurities etc., an incompressible ground state will be obtained. 

(b) The case v 2 = cooresponds to no electrons at all in the second layer, i.e. to the 
single layer case for which Jain proposed his ideas originally |l[]. For this case N 2 = and 
the second-layer coordinates Wi will be absent. Then ( |4.6|) , ( |3.6|) and (|3.11| ) reduce to 



4>v 1 (zi,z 2 , ,z Nl ) 



iVi 

n< 



Zi — Zj) 2si exp 



A'l 



4/ ; 



'pi 



-(2s 1 u 1 J2\zi\' 
i 

(zi,z 2 , ,z Nl ) 



(5.1) 



with V\ = , -jp - , This is just Jain's well known formula for the single layer spinless prob- 
lem. That our procedure for constructing non-unitary transformations to get flux-electron 
composites as developed in [17] will yield Jain's wavefunctions has also been pointed out by 



Wu and Yu [19] for the single layer case. Notice that the right hand side of eq.(5.1) contains 
a gaussian factor not included in Jain's version of this formula. Whether such a factor should 
be there or not just depends on the relative conditions under which the composite fermion 
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system is being compared to the electron system. The way Jain writes such an equation, the 
electron wavefunction -0 and the composite fermion wave function cj) correspond to the same 
magnetic field. Hence there is no relative gaussian factor between them. They also carry 
the same number of particles N\ . However they do correspond to different filling factors p\ 
and v\ and hence different densities from one another. This tacitly implies that in Jain's 
way of writing this relationship the two sides of the equation correspond to different areas 
(sample sizes). By contrast , we have defined the electron operator \l/ a snd the composite 
fermion operator Xa in the same domain, as is natural in a field theory. The sample areas 
are thus taken as equal. The total number of particles (and therefore the density) is also the 
same. The difference in filling factors is caused by the difference in effective magnetic fields, 
namely, B for the electron wavefunction i/j and B* (as given in eq (|3.10|) ) for the composite 
fermion state cj). Hence the gaussian factors (whose exponent is proportional to the magnetic 
field) will be different in ip and <p. The additional gaussian factor in ( |5.1|) compensates for 
this difference. To verify this note that eq( |3.10| ) reduces for the single layer case, to 



AB = B - B* = 1s x v x B (5.2) 

Recalling that B oc 4 we see that this difference AB corresponds precisely to the relative 
gaussian factor in (p.l|). 

(c) Returning to the double layer system, the various filling fraction possibilities con- 
tained in eq( |3.12D have been outlined at length by Lopez and Fradkin |2j . Of particular 



interest is the case of p\ = p 2 = lwith S\,S2 and n being arbitrary integers. Then the 



filling fractions (|3.12|) reduce to 

2s 2 + 1 - n 



(2si + 1) (2s 2 + 1) - n 2 
2si + 1 - n 
V2 ~ (2 fll + 1) (2s 2 + 1) - n" [ } 

Our theory then yields for the corresponding electronic wavefunction of this bilayer state 

with z/j, z/ 2 as given above, the formula (see |4.6|) : 

^ u u 2 {zi,z 2 , ,z Nl ; Wx,w 2 , ,w N2 ) 
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A'l 

IK 

i>j 



N 2 



JVi N 2 



/£>i Z n 



) 2si nK-^) 2s2 nn( 



Z r — We 



exp 



k>l 



Ni 



N 2 



^ i i 



W; 



h,i(zi,Z2, ,z Nl ; w 1 ,w 2 , ,w N2 ) (5.4) 



But 0i ; i is nothing but the wavefunction for unit filling factor in each layer, for each of which 
we can use the v — 1 Laughlin wavefunction, corresponding to an effective magnetic field of 
B* a = [ 1 — J2/3 Ra/3^/3 ) B = ( 1 — r a )B. This Laughlin wavefunction for unit filling in 
the first layer is 
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1> 



Laughlin 



]J(zi - Zj)exp 

i>3 



N 



47f ( ?^ 



(5.5) 



where 



1 _ eB\ (1-rQ 

Pi " 



(5.6) 



he I 2 

and similarly for the second layer in terms of the coordinates Wi . Inserting such a 0^1 into 
( ^4|) we get 



ij) UuV2 (z u z 2 , ,z Nl ; Wi,w 2 , ,w N2 ) 
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exp 
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I Ni N 2 

^tt 1 l 



UU 



(5.7) 



This is just the Halperin wavefunction for the ground state of the bilayer system || , 
with filling factors in the two layers given by (|5.3| ). Such a wavefunction was derived 



from a Chern-Simons field theoretic model by Ezawa and Iwasaki [16| earlier, but in a very 
different way. Their resuts are based on solutions to certain semiclassical self dual equations 
which require not only a mean-field approximation but also the neglect of e-e Coulomb 
interactions. They do treat interactions , but perturbativey, in the short distance limit. Our 
work here does use a mean field approximation in deriving results such as ( |5.4j ), but nowhere 
has the Coulomb term in the Hamiltonian dropped nor short distance appromation made. 
We also provide exact operator definitions of the composite particles. This wavefunction for 
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the case 2s\ + 1 = 2«2 + 2 (which is a special case of the above ) was also obtained by Lopez 
and Fradkin [0. They used their fermionic Chern-Simons theory [14]] suitably generalised 
to two-component wavefunctions and obtained the modulus of the wavefunction in the long 
wavelength limit. Our field theory, using the composite operators defined in ( |2.3| ) gives the 
full wavefunction including its modulus, phase and gaussian factors. No long wavelength 
approximation has been invoked by us. 

(d) Clearly our work is only a variant on the earlier work of Zhang et al, Read and Lopez 
and Fradkin cited above. The main advantage of our generalisation to non-unitary trans- 
formations and complex Chern Simons fields is that our composite operators incorporate all 
aspects of the Laughlin and Jain wavefunctions including the moduli \z{ — Zj\ of the corre- 
lations. These moduli contain the important zeroes which should be there in the presence 
of Coulomb repulsion, and also restore the wavefunctions to the lowest Landau level. The 
gaussian factors which should be present in the wavefunctions are also incorporated. Indeed, 
as our derivation shows these gaussian factors are essential for cancelling the imaginary part 
of our complex statistical gauge field in the mean field limit. True, some of these factors 
are obtained in the works of Zhang et al |P~3[ and Lopez and Fradkin M but only upon 



including fluctuations about the mean field. That they are present already at the mean field 
level in our operators indicates that our operators may be better candidates for flux-electron 
composite fields. 

(e) That we already get the full Laughlin and Jain wavefunctions at the mean field level 
raises hopes that corrections to these wavefunctions could be obtained even at the lowest 
order in fluctuations about the mean field. Unfortunately, this is where the non-unitary 
nature of our transformation (|2.3|) could create difficulties. Although the imaginary parts 
of our statistical field v(x) cancel in the mean field approximation they will be present 
away from mean field. Of course the full Hamiltonian ( |2.18| ) is hermitian as can be verified 
using ( |2.10| ), but its separation into a mean field part and a perturbation does not maintain 
hermiticity in each part. Standard perturbation techniques would have to be re-examined 
and modified to take this into account. These remarks hold not only for the present work 
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but also our earlier work with Sondhi [ HJ . For a discussion on how to go beyond mean field 
theory in such cases see Wu and Yu [|19j . 
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